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Abstract: We analyze nonlinear effects in optically driven arrays of 
nonlinear metallic nanoparticles. We demonstrate that such plasmonic 
systems are characterized by a bistable response, and they can support the 
propagation of dissipative switching waves (or plasmonic kinks) connecting 
the states with different polarization. We study numerically the properties 
of such plasmonic kinks which are characterized by a subwavelength extent 
and a tunable velocity. 
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1. Introduction 

When the size of conventional optical devices is reduced to nanoscale, the spatial confinement 
of light becomes inherently limited by diffraction. However, metal-dielectric structures allow 
to achieve a subwavelength light confinement with surface plasmon polaritons (SPPs), or plas- 
mons JT|[2][3). In particular, by using plasmons excited in an array of resonantly coupled metal 
nanoparticles JU [5], one can spatially confine and manipulate optical energy over distances 
much smaller than the wavelength. 

Because of the strong enhancement of the field induced by the excitation of SPPs and in- 
creased optical nonlinearity, surface plasmons can be employed for the realization of a va- 
riety of nonlinear optical effects. In particular, several nonlinear optical processes have been 
demonstrated in plasmonic nanostmctures, e.g., optical limiting and self-phase modulation in 
arrays of structured nanoparticles [0, second-harmonic generation and bistability in nanos- 
tructured metal films J7] |8] [9] and terahertz radiation generation in nonlinear metallic nan- 
odimers ifTOlfTTl . In addition, strong geometric confinement can boost efficiency of nonlinear 
optical effects, including the existence of subwavelength solitons in metal-dielectric multilay- 
ers lfl2l[T3lfT4l and arrays of metal nanowires lilSl . The solitons supported by such media result 
from a balance between tunneling of surface plasmon modes and nonlinear self-trapping. 

Here we study nonlinear effects in optically-driven arrays of nonlinear metal nanoparticles. 
In addition to the subwavelength plasmonic solitons analyzed earlier in other plasmonic sys- 
tems l!T2l n~3l IT4l IT31 . we demonstrate that an array of nonlinear metal nanoparticles can support 
the propagation of dissipative switching waves (or plasmonic kinks) connecting the states with 
different polarization. These kinks become possible due to a bistable nonlinear response of each 
nanoparticle 1 1 6 1 , and they have a subwavelength extent and a finite (or zero) velocity. 

The paper is organized as follows. In Sec. 2, we describe our nonlinear model and discuss its 
realistic parameters. Section 3 summarizes our results for the study of the basic properties of 
plasmonic kinks and their generation. Finally, Sec. 4 concludes the paper. 

2. Model and basic equations 

We consider a linear array of identical spherical metallic (silver) nanoparticles embedded into 
a SiC>2 host medium with permittivity £/,, and is exited by the optical field (see Fig. 1). 
We assume that the particle radii and the center-to-center distance are a = 10 nm and d = 
30 nm, respectively. Ratio a/d satisfies the condition ajd < 1/3, so that we can employ the 
point dipole approximation ifTTl . In the optical spectral range, a linear part of silver dielectric 
constant can be written in a generalized Drude form e^ g = £x, — (Op/ [a (to — z'v)], where = 
4.96, hto p — 9.54 eV, hv = 0.055 eV QU (hereinafter we accept exp(itot) time dependence); 
whereas dispersion of SiOi can be neglected since £/, ~ 2. 15 for wavelengths 350 - 450 nm |fl9l . 

Nonlinear dielectric constant of silver is = + X^l^n"" 1 ] 2 , where E„ is the local field 
inside n-th particle. We keep only cubic susceptibility due to spherical symmetry of particles. 



Fig. 1. Schematic sketch of an array of coupled metallic particles with a nonlinear response. 
Arrows indicate a profile of the polarizations for a typical kink configuration. 



According to the model suggested in Ref. 11201 and confirmed in experiment, 10 nm radii Ag 
spheres possess a remarkably high and purely real cubic susceptibility %^ ~ 3 x 10~ 9 esu, in 
comparing to which the cubic nonlinearity of SiC>2 is negligibly weak (~ 10~ 15 esu lED ). 

To study nonlinear dynamics of the particle polarizations, p„, we employ the model based 
on the dispersion relation method. This technique was firstly applied to investigation of light 
scattering by nonlinear metal nanodimers llTol ITTI and nonlinear effects in magnetic metama- 
terials Il22ll23l . and here we extend it to the case of an array of metal nanoparticle. According 
to this model, temporal dynamical response of metal nanoparticle arrays excited by the optical 
field with the frequency lying close to the frequency of the surface plasmon resonance of an 
individual particle, coq = (O p j\fe^ + 2e/„ can be described by the following system of coupled 
equations for slowly varying amplitudes of the particle dipole moments: 
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are in charge of dipole-dipole interactions between particles, P„ = 

p,i' ll VF 5 W 2 (£~ + 2£/ ! )£/ ! a 3 ) and = -3e h y^E^/yJ^ + 2 £/! ) 3 are di- 

mensionless slowly varying amplitudes of the particle dipole moments and external 
electric field, respectively, the indices "±" and 'if stand for the transverse and longitudinal 
components with respect to the array direction, 77 = - ^ (g ) 3 , |P„| 2 = l^l 2 + \P\\ 2 , 
7= v/(2<3Jo) + (^o fl ) 3 £ h I +2e/,) describes both thermal and radiation losses of particles, 
£0 = (Do/c-y/Eh, £1 — (ft) — a>o)/(0o and x = COot. Equations ([T} describe temporal nonlinear 
dynamics of an array of metal nanoparticles driven by arbitrary external optical field with the 
frequency (0 ~ cqq. We stress that model takes into account all particle interactions through the 
dipole fields, and it can be applied both to finite and infinite arrays. 




Fig. 2. (a) Dependence \Pq\ 2 on \Eq\ 2 at £1 = —0.1. Middle branch inside the bistability 
loop corresponds to the solutions unstable with respect to small perturbations. Arrows show 
the transitions of the system's state between the stable branches, (b) Profile of the function 
Re P£ for the stationary kink obtained by numerical simulations of Eq. QJ at £1 = —0.1 
and \E L \ 2 = 5.22 x 10~ 5 . 



3. Bistable response and plasmonic kinks 

In general, there are two fundamental nonlinear phenomena inherent in our system: modula- 
tional instability and bistability. Here we discuss the letter phenomenon. When the spatially 
distributed system demonstrates bistable response, it can support simultaneously two different 
stable states forming switching waves or kinks. However, to observe kinks and study its prop- 
erties, one should satisfy one more condition - the stability criterium, i.e. provide a range of the 
external field intensities (lying inside the bistability region) which is free from modulational 
instability. Our findings show that for an array of nanoparticles this condition is fulfilled for the 
transversal excitation only. Therefore, we, next, consider just this case. 

To begin with, we excite the infinite array by an homogenous electric field with the transver- 
sal polarization: E„ = (Eq,0). In this case, all particle dipole moments remain the same, 
P„ = Pq, and the stationary solution of Eq. (fTJ can be written as follows 

(-iy+Cl + pAj + \I$A P£ = Efr, (2) 

where 

/ 1 ik d (k d) 2 \ 
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Transition from G,j„, to Aj has been made via the replacement \n — m\ = j and taking into 
account symmetry structure of the series. When 



Q. < -Re £ Aj - V3 ^y-Im £ t Ajj , 



the polarization Pq becomes a three-valued function of Eq leading to bistability (see Fig.[2|a)). 

To analyze switching waves, we perform numerical simulations of Eq. (QJ for a finite array 
(with 200 nanoparticles) at initial conditions corresponding to different branches of the homo- 
geneous stationary solution (01 over different parts of the array. In Fig. |2b) the profile of the 
stationary kink is presented. However, our calculations show that kinks can move as well. We 
have investigated impact of the field intensity on temporal kink dynamics. Characteristic results 
are summarized in Figs.[3la,c). To simplify analysis, we fixed Q. and set E^ as a step-like func- 
tion of time, as shown in Figs.[5Jb,d). Constant Eq results in the constant kink velocity, v. That 
is why a step-like time dependency of Eq allows us to control v overall bistability zone. 
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Fig. 3. (a) (Media 1) and (c) (Media 2) demonstrate dynamics of Re P t j obtained numeri- 
cally from Eq. |QJ at £2 = —0.1. (b,d) External light intensities corresponding to the cases 
(a) and (c), respectively. Dashed lines mark the bistability region. 
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Fig. 4. Normalized kink velocity vs. the intensity of applied field \Eq\ at fl = —0.1. 



Figure |3ja) demonstrates dynamical behavior of the system when Eq increases from the 
lower to the upper threshold of the bistability region, as shown in Fig. Ob). In this case the 
kink state changes from standing to moving, and for \Eq\ 2 lying close to the middle of the 
bistability region the second kink is appeared at the left edge of the array, moving towards the 
first one. This is the straight consequence of cooperation between edge effects and modulational 
instability, which, in the end, crushes both kinks. 

To avoid appearance of modulational instability, we also carry out the simulation when Eq 
decreases from the middle to the lower threshold of the bistability region. The result is presented 
in Fig. EJc). Here the kink velocity changes from the finite value to zero, and modulational 
instability growth is not observed. 

Figure|4]shows the kink velocity as a function of the external field intensity \Eq | 2 . We notice 
that there exists a finite range of the field intensities where the kink does not move, i.e. u = 0. 
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Fig. 5. Generation of a plasmon soliton via interaction of two kinks with the opposite 
polarities, (a) Initial condition used for numerical simulations of Eq. ([TJ at l^l 2 = 9.27 x 
10~ 5 and £2 = —0.1. (b) (Media 3) The frozen frame of Re showing the profile of the 
subwavelength dissipative plasmon soliton formed by two interacting kinks. 



Generally speaking, this is a common property of nonlinear discrete systems |24| where the 
effective periodic potential created by the lattice requires a finite value of the applied external 
force to start kink's motion. The existence of the kinks with zero velocity suggest a straight- 
forward way to construct spatially localized modes by taking two standing kink of the opposite 
polarity. The resulting localized modes form the continuous families, and they resemble the 
similar localized modes discussed earlier for metal-dielectric multilayers lfT2] [PJl [141 and ar- 
rays of metal nanowires |fT31 . On the other hand, we may expect generating such subwavelength 
plasmon solitons by considering interaction of two slowly moving kinks of the opposite polar- 
ity. 

Here, we demonstrate this latter approach. We study numerically the interaction of two kinks 
of the opposite polarities, and select the value of the external field to realize the kinks with 
the near-zero velocity v [see Fig. HJa)]. We find that the kinks start moving slowly towards 
each other, then stop and form a stable dissipative plasmon soliton, as shown in Fig. [3b). The 
generated plasmonic mode is localized only on 6 particles, and its width is about 150 nm. 
Because £2 = —0.1 corresponds to the radiation wavelength A = 440 nm, this plasmon soliton 
is truly subwavelength because it extends for 0.34A. 

It remains open if the kinks motion may demonstrate bistability, as was recently revealed for 
a similar nonlinear system with local coupling [23]. This will be a subject for future studies. 

4. Conclusions 

We have demonstrated that an array of interacting nonlinear metal nanoparticles can support 
a novel type of nonlinear localized modes, plasmonic kinks, which describe switching waves 
connecting two different states of polarization of the metal particles. We have generated various 
types of such kinks, analyzed numerically their velocity as a function of the external field 
intensity, and revealed that the array of metal nanoparticles can support deeply subwavelength 
dissipative plasmon solitons which in our case can be regarded as a bound state of two kinks of 
the opposite polarity. 
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